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ABSTRACT 

This  report  describes  an  implementation  of  a variation  of  the 
hidden-line  algorithm  developed  by  P.  Loutrel.  The  new  algorithm 
produces  hidden-line  drawings  of  scenes  composed  of  opaque  polyhedra 
as  well  as  of  some  more  general  types  of  flat-faced  objects  not 
covered  by  the  original  Loutrel  algorithm.  The  program  is  written  in 
Fortran  IV  and  has  been  implemented  in  a batch  version  on  a CDC  6600 
computer  and  in  an  interactive  version  on  an  .ADAGE  AGT-30  interactive 
graphics  computer. 
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1.  INTRODUCTION 


Given  a computer  description  of  a three-dimensional  scene  and  the 
coordinates  of  a vantage  point , a hidden-line  algorithm  determines  the 
parts  of  the  scene  which  are  visible  from  the  vantage  point. 

In  this  report  an  implementation  of  such  an  algorithm  is  described. 
The  algorithm  is  that  originally  developed  by  P.  Loutrel  and  described 
in  a journal  article*  as  well  as  in  his  doctoral  thesis**.  The  original 
algorithm  accepted  only  true  polyhedral  objects,  either  convex  or  concave 
and  either  simply  or  multiply  connected.  The  algorithm  implemented  here 
extends  the  set  of  acceptable  objects  to  all  objects  bounded  by  planar 
surfaces.  This  includes  objects  whose  faces  may  contain  one  or  more 
holes  (either  intruding  or  extruding).  However,  the  extended  algorithm 
cannot  handle  non-closed  objects  or  objects  bounded  by  quadric  or  other 
non-planar  surfaces. 

The  implemented  program  called  PD  RAW  is  written  in  FORTRAN  IV.  There 
are  two  versions  of  the  program:  batch  and  interactive. 


* Loutrel,  ?.  P. , "A  Solution  to  the  Hidden-Line  Problem  for  Computer- 
Drawn  Polyhedra,"  I5EE  Transactions  or.  Ccmruters,  Voi.  C-i  No.  ?, 
March  1970,  pp.  205-213. 

**  Loutrel,  P.  P. , "A  Solution  to  the  Hidden-Line  Problem  for  Computer- 
Drawn  Polyhedra,"  doctoral  dissertation.  Department  of  Electrical 
Engineering,  New  York  University,  September  196?.  (NTI3  accession 
number:  No8  13  S30). 
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The  batch  version  was  implemented  on  a CDC  66 00  computer, 
version  of  the  program  reads  all  commands,  object  descriptions  and 
vantage  point  specifications  from  an  input  data  file  and  produces 
drawings  on  a digital  plotter. 

The  interactive  version  was  implemented  on  an  ADAGE  AGT-30  inter- 
active graphics  computer.  This  version  of  the  program  reads  object 
descriptions  from  data  files  but  program  control  and  vantage  points  are 
specified  by  interactive  graphics  I/O  devices.  The  drawings  produced 
by  the  program  are  made  on  a CRT  screen  and  hard  copies  can  be  obtained 
on  a dot  matrix  plotter  or  an  incremental  plotter. 

During  the  last  several  years  a number  of  other  hidden-line  or 
hidden-surface  algorithms  have  been  developed.  An  extensive  survey  of 
these  algorithms  is  given  in  [3]. 
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2.3  An  Extension  of  the  Algorithm 

The  original  algorithm  was  extended  to  accept  objects  that  have  one 
or  more  face-interior  holes  [Fig.  2.2]  or  protrusion  [Fig.  2.3].  When 
an  object  has  a face-interior  hole  or  protrusion,  the  face  where  the  hole 
or  protrusion  occurs  will  be  multiply  connected;  that  is,  it  will  have 
a hole  in  it.  We  shall  refer  to  a hole  in  a face  as  an  "antiface".  When- 
ever a vertex  is  behind  a visible  face,  the  order  of  invisibility  is 
increased  by  one.  However,  when  a vertex  lies  behind  an  antiface,  the 
order  of  invisibility  is  decreased  by  one.  This  is  done  because  each 
antiface  indicating  the  presence  of  a face-interior  hole  or  protrusion 
is  inside  a larger  true  face*.  For  the  face  the  order  of  invisibility 
is  increased;  for  the  antiface  the  order  of  invisibility  is  decreased*. 
Therefore,  the  order  of  invisibility  of  a vertex  inside  an  antiface  does 
not  change  after  both  the  face  and  the  antiface  have  been  processed. 


* An  antiface  is  thus  bounded  by  a chain  of  edges  that  is  di.sj  :int 
from  the  chain  bounding  the  containing  face. 


3.  THE  PD RAW  ROUTINES 


In  this  section  all  routines  in  PDRAW  are  described.  PDRAW  con- 
sists of  a main  program,  9 subroutines  and  h functions.  The  chart 
of  Fig.  3.1  gives  their  connections. 


3.1  Program  PDRAW 

PDRAW  is  the  main  program  which  controls  the  execution  of  all 
routines.  It  opens  the  plotting  output,  moves  the  pen  to  the  first 
drawing,  moves  the  pen  between  drawings  and  closes  the  plotting  output. 
PDRAW  also  reads  from  a data  file  all  the  data  delimiters  and  vantage 
point  specifications,  analytes  them  and  calls  the  main  functional  sub- 
routines . 


3 . 2  Subroutine  SETUP 

SETUP  is  a subroutine  which  initializes  the  data  structure  of  the 
program.  Pointers  to  the  various  sections  of  the  data  structure  are  se‘ 
to  point  to  the  first  available  locations,  and  the  number  of  objects  is 
set  to  zero.  The  data  structure  of  the  program  is  described  in  Section 


3 . 3  Subroutine  INPUTM 

This  subroutine  reals  the  description  of  an  object  from  the  input 
data  file  to  the  data  structure  of  the  program.  The  format  of  an  objec' 
description  is  described  in  Section  ;.  Whenever  an  attempt  is  made  to 
read  more  data  into  the  data  structure  than  the  allowable  maximum,  an 
error  message  is  written,  the  plotting  output  is  closed  and  program 


execution  is  stopped. 


3 . U Subroutine  PRSCCM 

This  subroutine  performs  the  preliminary,  vantage-point-independent 
classification  of  each  object.  It  is  called  once  for  each  object 
immediately  following  a call  to  INPUTM . PRECOM  computes  the  following: 

a)  an  outward  normal  vector  to  each  face  of  the  object 
(the  method  is  described  in  Section  2.2) 

b)  classification  of  each  edge  as  either  type  H2  or  H3 
[l,  page  208] 

c)  classification  of  the  object  as  either  convex  or  concave. 

3 . 5 Subroutine  TRANS 

This  subroutine  transforms  all  vertices  of  all  objects  from  the 
object  coordinates  to  the  vantage  point  coordinates  and  then  to  the 
picture  plane  coordinates.  TRANS  uses  either  perspective  or  orthographic 
projection.  All  picture  plane  coordinates  are  scaled  down  to  a 
(-1.0,  1.0)  square  for  further  computations.  The  transformation  process 
is  described  in  Appendix  A. 


3.0  Subroutine  REMOVE 


■ controlling 

the  actual  tas 

k of 

hidden- 

■isibility  oi 

all  faces  by  c 

alls 
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vertices  are 

classified  by 

call 

ing  sub- 

routine  CLASS.  Visible  edges  in  a non-intersecting  convex  polyheir- n* 
are  determined  in  this  subroutine.  Each  scene  is  diviiei  into  sets  of 
potentially  intersecting  object#  and  passed  to  subroutine  LOCATE  f?r 
processing. 


* Intersections  are  determined  for  projecti  ons  ejects  in*  * : 
picture  plane. 
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3.7  Function  IVFACE 

'This  is  an  integer  function  which  returns  the  value  of  +1  if  a 
face  is  visible  from  a given  point  and  the  value  of  -1  if  the  face  is 
invisible  from  a given  point.  IVFACE  is  invoked  at  two  points  in  the 
program.  First,  it  is  used  to  determine  whether  each  face  is  invisible 
or  potentially  visible  from  a vantage  point.  Second,  in  order  to  compute 
the  order  of  invisibility  of  a vertex,  this  function  is  used  to  determine 
whether  a face  is  visible  from  a vertex,  which  is  identical  to  deter- 
mining whether  the  face  is  in  front  of  or  behind  the  vertex.  ’This 
second  use  of  the  function,  however,  uses  only  faces  which  were  deter- 
mined to  be  visible  in  the  first  use  of  the  function,  since  invisible 
faces  do  not  change  the  order  of  invisibility  of  a vertex. 

3.8  Subroutine  MRANGS  : 

This  subroutine  computes  an  enclosing  envelope  for  each  object  in 

the  current  scene.  This  is  done  by  searching  for  the  maximum  and  minimum 
x and  y picture  coordinates  of  each  object. 

3 . 9 Subroutine  FRAilGE 

This  subroutine  computes  an  enclosing  envelope  for  each  visible 
face  of  each  object.  This  is  done  by  searching  for  the  maximum  and  mini- 
mum x and  y picture  plane  coordinates  of  each  visible  face  in  a scene. 

Although  the  enclosing  envelopes  of  visible  faces  are  not  required  by 
the  basic  algorithm,  they  are  used  to  speed  up  substantially  the  compu- 
tations of  the  order  of  invisibility  of  vertices  as  well  as  the  searching 
for  edge  intersections  performed  in  subroutines  COLECT  and  IORDER. 
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3.10  Subroutine  CLASS 

This  subroutine  performs  the  vantage-point-dependent  classifica- 
tions. Each  edge  in  a scene  is  classified  as  either  type  HI,  H2,  H3 
or  HU  [l,  pp.  208-209].  Each  vertex  is  classified  as  either  a boundary 
or  an  inner  vertex  [l,  pp.  210-211].  This  subroutine  is  called  each 
time  a set  of  potentially  intersecting  objects  is  found.  If  it  is 
called  with  only  a single  object  which  is  convex,  the  vertex  classifi- 
cation is  omitted.  This  subroutine  identifies  all  invisible  edges 
(type  HI  and  H2). 


3.11  Function  IORDER 

This  function  returns  as  its  integer  value  the  order  of  invisibility 
[l,  page  209]  of  a vertex  passed  to  it.  The  order  of  invisibility  is  a 
non-negative  integer  corresponding  to  the  number  of  potentially  visible 
faces  which  are  in  front  of  the  vertex,  i.e.  between  the  vantage  point 
and  the  vertex.  Note  that  invisible  faces  are  not  counted. 

3.12  Subroutine  LOCATE 

REMOVE  passes  to  LOCATE  a set  of  objects  in  a scene  which  potentially 
intersect  in  the  picture  plane  [l,  page  211 ].  This  routine  processes 
in  two  passes  all  H3  and  Hi  type  edges.  In  the  first  pass  all  K3  edges 
are  processed;  in  the  second  pass  all  Hi  edges  are  processei.  Each  edge 
that  is  processed  has  its  initial  and  final  coordinates  put  into  a stack. 
The  order  of  invisibility  of  the  initial  vertex  of  the  edge  is  put  into 
another  stack.  These  stacks  are  then  passed  to  subroutine  CCLECT  which 
computes  any  possible  intersection  of  this  edge  with  any  H3  or  Hi  type 
edge  in  the  set  of  objects  passed  to  LOCATE. 
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3.13  Subroutine  COLECT 

This  subroutine  computes  all  possible  intersections  that  one  edge 
may  have  with  another  [l,  pp.  209-210].  The  initial  and  final  coordi- 
nates of  an  edge  are  passed  to  COLECT  from  LOCATE.  COLECT  uses 
functions  ITEST  and  INTERS  to  find  any  possible  intersections  and  adds 
the  coordinates  of  any  found  intersection  into  a stack.  It  also  com- 
putes whether  the  processed  edge  is  coming  from  behind  a face  or  is 
going  behind  a face  at  the  point  of  the  intersection.  This  change  in 
the  order  of  invisibility  along  the  edge,  which  is  either  +1  or  -1,  is 
also  put  into  a stack. 

3- Ik  Function  ITEST 

ITEST  is  an  integer  function  which  tests  two  edges  for  a valid 
intersection.  It  returns  the  value  of  +1  if  a valid  intersection  is 
found;  otherwise  the  value  of  0 is  returned.  ITEST  is  also  used  to  find 
an  intersection  of  a line  going  from  a vertex  to  infinity  with  an  edge 
during  computations  of  the  order  of  invisibility  of  a vertex. 

3-15  Function  INTERS 

INTERS  is  called  from  COLECT  after  ITEST  found  that  two  edges  have 
a valid  intersection  and  computed  the  x and  y coordinates  of  the  inter- 
section. INTERS  determines  whether  at  the  point  of  the  intersection , the 
processed  edge  is  in  front  of  the  other  edge  or  not  [l,  page  210’.  If 
it  is  behind,  it  means  that  the  processed  edge  is  either  coming  *'rcn 
behind  a face  or  is  going  behind  a face  and  INTERS  returns  the  value  of 


+1 , otherwise  a zero  is  returned. 
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3.16  Subroutine  EMIT 

This  subroutine  is  passed  a stack  containing  the  initial  and  final 
coordinates  of  an  edge  as  well  as  the  coordinates  of  all  intersections 
this  edge  has  with  all  other  edges.  A stack  containing  the  initial 
order  of  invisibility  as  well  as  the  changes  in  the  order  of  invisibi- 
lity at  each  intersection,  either  +1  or  -1,  is  also  passed  to  EMIT. 

EMIT  sorts  out  the  intersections  by  their  increasing  distance  from  the 
initial  vertex  and  then  causes  the  segments  of  the  edge  to  be  drawn  as 
either  visible  or  invisible  vectors. 
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4.  THE  DATA  STRUCTURE 

In  this  section  the  data  structure  of  the  program  is  described. 
All  the  data  describing  objects  as  -well  as  several  parameters  used 
throughout  the  program  are  stored  in  common  blocks.  The  portions  of 
the  data  structure  which  are  vantage-point  independent  are  marked 
with  asterisks.  A chart  of  the  data  structure  is  given  in  Fig.  4.1. 

4 . 1 Common  Block  MODEL 

This  common  block  contains  the  data  describing  the  objects  that 
have  been  entered  into  the  program.  All  arrays  in  this  block  are 
dimensioned  to  MAXM  which  is  the  maximum  number  of  objects  in  a scene 
that  can  be  entered  into  the  program. 

MODEL ( MAXM ) * contains  pointers  to  the  first  face  of 
an  object.  The  pointer  of  each  object  is  set  by 
subroutine  INPUT  just  before  the  object  is  read  in. 

MTYPE(MAXM)*  contains  the  type  of  an  object: 

convex  object  0 

concave  object  1 

The  type  is  set  by  subroutine  PRECCM  after  the 
preliminary  classification  of  all  edges  as  either 
type  H2  or  H3.  If  there  are  no  edges  of  type  H2, 


■■■■ 


- 1 6 - 

IFRSTV(MAXM)*  contains  pointers  to  the  first  vertex  of 
of  each  object.  The  pointer  of  an  object  is  set  by 
subroutine  INPUTM  just  before  the  object  is  read  in. 

ILASTV(MAXM)*  contains  pointers  to  the  last  vertex  of 
each  object.  The  pointer  of  an  object  is  set  by 
subroutine  INPUTM  after  the  last  vertex  of  the  object 
was  read  in. 

XMAXM ( MAXM ) , XMINM ( MAXM ) , YMAXM ( MAXM ) and  YMINM ( MAXM ) 

contain  the  maximum  and  minimum  x and  y picture  plane 
coordinates  of  each  object.  These  values  are  computed 
for  all  objects  in  subroutine  MRANGE. 

MLIST(MAXM)  is  a processing  status  list  which  contains 
pointers  to  a set  of  objects  that  have  potential 
intersections  in  the  picture  plane.  These  pointers 
are  set  in  subroutine  REMOVE. 

MSTATE(MAXM)  contains  the  processing  status  of  each  object 


in  a scene: 

not  processed  0 

add  to  the  processing  list  ...  1 

in  the  processing  list  2 

processed  3 
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4.2  Common  Brock  FACE 

This  common  block  contains  the  data  describing  the  faces  of  objects 
entered  into  the  program.  All  arrays  in  this  block  are  dimensioned  to 
MAXF,  which  is  the  maximum  number  of  faces  that  can  be  entered  into  the 
program . 

IFACE(MAXF)*  contains  pointers  to  the  first  edge  of  a 

face.  The  pointer  of  each  face  is  set  by  subroutine 
INPUTM  just  before  a face  is  read  in. 

FNV(MAXF,3)*  contains  the  three  components  of  an  out- 
ward normal  vector  of  each  face.  The  vector  is 
computed  in  subroutine  PRECOM. 

IVPNT ( MAXF ) * contains  pointers  to  the  firsu  vertex  of 
each  face.  The  pointer  is  set  by  subroutine  PRECOM. 

NEXTF(MAXF)*  contains  pointers  to  the  next  face  of  an 
object.  All  pointers  are  initialised  by  subroutine 
SETUP.  End-of-list  pointers  are  inserted  by  sub- 
routine INPUTM  after  the  last  face  of  an  object  was 
read  in. 

IrTYPE ( MAXF ) * contains  the  type  of  each  face; 

inner  face  description  0 

outer  face  description  1 

The  type  of  each  face  is  read  from  the  input  data 


file  by  subroutine  INPUTT!. 
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ICF(MAXF)  indicates  the  visibility  of  each  face: 

invisible  face  -1 

potentially  visible  face  +1 

The  visibility  is  computed  in  function  IVFACE  and 
stored  to  ICF  by  subroutine  REMOVE. 

XMAXF(MAXF),  XMINF(MAXF),  YMAXF ( MAXF ) and  YMINF(MAXF) 

contain  the  maximum  and  minimum  x and  y picture  plane 
coordinates  of  each  potentially  visible  face.  These 
values  are  computed  for  all  potentially  visible 
faces  by  subroutine  FRANGE. 

ZFNV(MAXF)  is  used  in  computing  the  visibility  of  faces 
when  the  orthographic  projection  is  selected.  It 
contains  the  dot  product  of  the  o row  of  the  trans- 
formation matrix  with  a face  normal  vector.  A 
negative  value  indicates  an  invisible  face;  a non- 
negative face  value  indicates  a potentially  visible 
face . 

h . 3 Common  Block  EIOE 

This  common  block  contains  the  data  describing  the  edges  of  objects 
entered  into  the  program.  All  arrays  in  this  block  are  dimensioned  to 
MAXE  which  allows  MAXE/2  edges  to  be  entered  into  the  program. 

IEDGE(MAXE,k)*  contains  the  following  four  pointers 
for  each  edge: 

1.  Pointer  to  the  initial  vertex  of  the  edge. 

final  vertex  of  the  edge. 


2.  Pointer  to  the 
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3.  Pointer  to  the  face  on  the  right  of  the  edge. 

L.  Pointer  to  the  face  on  the  left  of  the  edge. 

Pointers  1 and  2 are  read  in  by  subroutine  INPUTM  from 

the  input  data  file.  Pointers  3 and  4 are  set  by  subroutine 
INPUTM  after  all  the  data  was  read  it . 

IHTYPE(MAXE)*  contains  the  vantage-point-independent  classifi- 
cation of  each  edge  as  either  type  H2  or  type  H3.  This 
classification  is  done  in  subroutine  PRECOM. 

NEXTE(MAXE)*  contains  pointers  to  the  next  edge  in  a face  of 
an  object.  All  pointers  are  initialized  by  subroutine 
SETUP.  End-of-list  pointers  are  inserted  by  subroutine 
INPUTM  after  a face  was  read  in. 

IH(MAXE)  contains  vantage-point-dependent  classification 
of  each  edge  as  type  HI,  H2,  H3  or  HU.  This  classi- 
fication is  done  by  subroutine  CLASS. 

U . U Common  Block  VERTEX 

This  common  block  contains  the  data  describing  the  vertices  of  objec 
entered  into  the  program.  All  arrays  in  this  block  are  dimensioned  t,o 
MA XV  which  allows  MAXV  - 1 vertices  to  be  entered  into  the  program.  The 
last  location  is  used  by  the  program  for  a point  at  "infinity". 

VERTEX ( MAXV, 3)*  contains  the  coordinates  of  vertices 
described  in  the  object  coordinate  system.  The 
coordinates  are  read  in  from  the  input  data  file 


by  subroutine  INPUTM. 
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VTRANS(MAXV,3)  contains  the  coordinates  of  vertices  in  the 
vantage  point  coordinate  system.  The  transformation 
from  the  object  coordinate  system  to  the  vantage  point 
coordinate  system  is  done  by  subroutine  TRANS. 

V2D(MAXV,2)  contains  the  coordinates  of  vertices  in  the 
picture  plane  coordinate  system.  The  transformation 
from  the  vantage  point  coordinate  system  to  the  picture 
plane  coordinate  system  is  done  by  subroutine  TRANS. 

XOV(MAXV)  cont  ains  the  order  of  invisibility  of  a vertex. 

It  is  computed  by  function  IORDER. 

IVTYPE(.MAXV)  contains  the  classification  of  a vertex  as 
either  a boundary  vertex  (non-zero  value)  or  an 
inner  vertex  (zero  value). 

U . 5 Common  Block  STACK 

This  common  block  contains  information  about  intersections  that  an 
edge  may  have  with  any  other  edge.  All  arrays  in  this  block  are  dimen- 
sioned to  MANS.  This  allows  each  edge  to  have  up  to  MAXS  - 2 inter- 
sections, since  two  locations  are  used  for  the  initial  and  final  vertices 
of  the  edge. 

I3TPNT  is  a pointer  to  the  top  of  the  stack. 

STACK(.MAX3,  2)  contains  the  picture  plane  x and  y coordi- 
nates of  all  edge  intersections.  Locations  1 and 
2 contain  the  picture  plane  coordinates  of  the 
initial  and  final  vertices  of  "he  edge. 


L. 
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KS(MAXS)  contains  the  increment  of  the  order  of  invisi- 
bility of  each  intersection.  This  value  may  be  either 
+1  or  -1  indicating  that  the  edge  is  going  behind  a 
face  or  coming  from  behind  a face.  Locations  1 and 
2 of  the  stack  contain  the  order  of  invisibility  of 
the  initial  and  final  vertices  of  the  edge. 

DIST(.MAXS)  is  used  for  sorting  out  edge  intersections  by 
their  increasing  distance  from  the  initial  vertex 
of  the  edge.  The  sorting  is  done  in  subroutine 
EMIT  before  the  subroutine  emits  segments  of  the 
edge  as  either  visible  or  invisible  vectors. 

h.6  Common  Slock  FREE 

This  common  block  contains  pointers  to  the  next  available  locations 
in  the  blocks  MODEL,  FACE,  EDGE  and  VERTEX. 

IMFREE*  is  a pointer  to  the  next  available  location  in 
block  MODEL.  It  is  initialized  to  1 by  subroutine 
SETUP  and  incremented  by  1 in  subroutine  IIIPUTM 
before  a new  object  is  read  in. 

IEEREE*  is  a pointer  to  the  next  available  location  in 
block  FACE.  It  is  initialized  to  1 by  subroutine 
SETUP  and  incremented  by  1 in  subroutine  INPUTM 
before  a new  face  is  read  in. 

IEFP.EE*  is  a pointer  to  the  next  available  location  ir. 
block  EDGE.  It  is  initialized  to  1 by  subroutine 
SETUP  and  incremented  by  1 in  subroutine  INPUTM 
before  a new  edge  is  stored  in. 
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IVFREE*  is  a pointer  to  the  next  available  location  in 
block  VERTEX.  It  is  initialized  to  1 by  subroutine 

I 

SETUP  and  incremented  by  1 in  subroutine  INPUTM 
before  a new  vertex  is  read  in. 

NMODEL*  is  the  number  of  objects  in  the  scene  in  the 
program.  It  is  set  to  0 by  subroutine  SETUP  and 
incremented  by  1 in  subroutine  INPUTM  before  a 
new  object  is  read  in. 

NMSET  is  the  number  of  objects  in  the  scene  which  have 
potential  intersections  in  - he  picture  plane  and 
which  are  currently  in  the  list  MLIST.  The  number 
is  computed  in  subroutine  REMOVE. 

U . 7 Common  31ock  PARAMS 

This  common  block  contains  several  parameters  which  are  used  t 
out  the  program. 

VP(3)  contains  the  current  vantage  point  x,  y and  z 
coordinates  given  in  the  object  coordinate  system. 

They  are  read  from  a vantage  point  specification 
card  by  the  main  program  PDRAW.  In  seme  routines 
these  variables  are  set  equivalent  to  X0,  Y0  and  20. 

ORTKO  is  a logical  variable  which  is  set  to  .TRUE,  in 
the  main  program  PDRAW  when  the  orthographic 
projection  is  selected. 

ZR(3)  contains  the  z row  of  the  transformation  matrix, 
computed  in  subroutine  TRANS. 


L. 
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DD  contains  the  transformation  value  D,  the  distance  of 
the  vantage  point  from  the  origin  of  the  object 
coordinate  system,  computed  in  subroutine  TRANS. 

ZMAX  contains  the  picture  scale  of  the  current  image. 

It  is  computed  in  subroutine  TRANS . 

RADIUS  is  the  maximum  distance  of  a vertex  from  the 
origin  of  the  object  coordinate  system.  It  is 
computed  in  subroutine  INPUTM. 

EPS1,  EPS2,  EPS3  are  small  numbers  which  are  used  as 
tolerances  in  some  computations.  All  must  be 
positive  numbers. 

IENPNT  is  a pointer  to  the  currently  processed  edge. 

N1  is  a pointer  to  the  initial  vertex  of  the  currently 
processed  edge.  It  is  set  in  subroutine  LOCATE 
and  used  in  programs  COLECT,  ITEST,  and  INTERS. 

N2  is  a pointer  to  the  final  vertex  of  the  currently 
processed  edge. 

ISIPNT  is  a pointer  to  an  e dge  which  may  have  an  inter- 
section with  edge  IENPNT  and  is  currently  tested. 

11  is  a pointer  to  the  initial  vertex  of  edge  IEIPNT. 

It  is  set  in  subroutine  COLECT  and  used  in  programs 
ITEST  and  INTERS. 

12  is  a pointer  to  the  final  vertex  of  the  edge  IEIPNT. 
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X is  the  x coordinate  in  the  picture  plane  of  the 
intersection  (if  it  exists)  of  edges  IENPNT  and 
IEIPNT.  It  is  computed  in  function  ITEST  and 
used  in  programs  COLECT  and  INTERS. 

Y is  the  y coordinate  of  the  intersection. 

DASHED  is  a logical  variable  set  to  .TRUE,  by  PDRAW 
whenever  the  hidden  edges  are  to  be  drawn  as 
dashed  lines. 

SCALE  sets  the  plot  output  for  a maximum  plotting  area 
of  2*  SCALE  by  2*  SCALE  inches. 

L.8  Common  31ock  DIMENS 

This  common  block  contains  the  dimensions  of  the  arrays  used  in  the 
data  structure. 


MAXM* 

contains 

the 

dimension 

of 

arrays 

in 

block 

MODEL. 

MAX?* 

contains 

the 

dimension 

of 

arrays 

in 

block 

PACE. 

MAXE* 

contains 

the 

dimension 

of 

arrays 

in 

block 

EDGE. 

MAXV* 

contains 

the 

dimension 

of 

arrays 

in 

block 

VERTEX 

MAXS* 

contains 

the 

dimension 

of 

arrays 

in 

block 

bTAC K . 

5 . DATA  FORMAT 


In  this  section  the  format  of  data  files  read  by  the  program  is 
described.  A data  file  contains  three  types  of  cards  - or  card  images. 

delimiter  cards  - delimit  sections  of  a data  file. 

object  cards  - contain  object  description. 

vantage-point  cards  - contain  vantage  point  description. 

5 . 1 The  Structure  of  Data  Files 

Input  data  to  PDRAW  consists  of  a description  of  one  or  more  scenes. 
Each  scene  contains  one  or  more  objects.  The  last  object  in  a scene  is 
followed  by  one  or  more  vantage  point  specifications.  The  format  of 
description  of  an  object  is  given  in  5-2.  A scene,  an  object  and  a 
group  of  vantage  points  are  delimited  by  special  cards  - delimiter  cards  - 
which  contain  a digit  in  column  one  as  follows: 

1 Scene  header;  indicates  that  a scene  follows;  columns 
2-1+2  may  contain  the  name  of  the  scene  or  any  comment 
which  is  printed. 

2 Object  header;  indicates  that  a description  of  an  object 
follows  (in  the  format  described  in  5.2);  columns  2-1+2  may 
contain  the  name  of  the  object  or  any  comment  which  is 
printed. 

6 Vantage-point  header;  indicates  that  vantage  point 

cards  follow;  this  card  is  placed  following  the  descrip- 
tion of  the  last  object  in  a scene. 

9 End-of-scene  card.  This  card  follows  the  last  vantage 
point  card  and  indicates  the  end  of  the  current  scene. 
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It  may  be  followed  by  an  end-of-file  mark  or  by  the 
header  of  the  next  scene.  The  last  scene  in  a file  need 
not  be  followed  by  an  end-of-scene  card. 

The  sice  of  the  input  data,  i.e.  the  number  of  objects,  faces,  edges  and 
vertices  in  a scene  is  limited  by  the  size  of  the  data  structure.  An 
error  message  is  printed  when  the  permitted  size  of  the  data  structure 
is  exceeded,  indicating  the  part  of  the  data  structure  that  was  exceeded. 
The  program  does  not  limit  the  number  of  scenes  or  vantage  points  in  a 
data  file. 


5.2  Format  of  Object  Specifications 

In  order  to  enter  a planar-faced,  solid  object  into  the  program,  a 
description  of  the  object  must  be  created.  Such  a description  consists 
of  a list  of  the  vertices,  followed  by  a list  of  the  faces  of  the  object. 
5.2.1  Vertices 

Each  vertex  V.  is  defined  bv  three  coordinates 
i 

V.  (x.  ,y  ,z  ) given  in  the  object  coordinate  system.  The  index  ' i ' may 
be  arbitrarily  assigned  to  a vertex.  It  is,  however,  assumed  that 
i = 1,2,3, M where  M is  the  total  number  . f verti  res.  f ra&t 


the  list 


M 


x v z 
1 ^l  1 


X2  V2  Z2 


13  (total  number  of  vertices 

3F8.2  ^ first  vertex 

3F8.2  (second  vertex 


7M  ZM 


3Fa . 2 


(last  vertex 
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5.2.2  Faces 

Each  face  of  a planar-faced,  solid  object  is  a polygon. 

The  program  allows  faces  to  have  one  or  more  holes  in  them,  that  is,  to 
contain  antifaces.  Therefore,  to  describe  a face,  a list  of  vertices 
defining  the  outer  edges  as  well  as  lists  defining  any  possible  holes 
or  protrusions  must  be  given.  Each  such  list  contains  the  indices  ' i ' 
of  vertices  as  defined  in  5.2.1.  The  direction  in  which  vertices  are 
listed  for  a face  is  such  that  the  interior  lies  always  to  the  right 
relative  to  an  observer  outside  of  the  object.  For  an  antiface,  the 
direction  is  taken  so  that  the  interior  lies  to  the  left.  Each  list  of 
vertices  must  be  closed;  i.e.,  the  first  vertex  must  appear  also  as  the 
last  one.  Note  that  all  the  vertices  and  edges  of  a face  (including 
those  of  all  its  antifaces,  if  any)  must  lie  in  a plane. 


There 

are 

N lists  of  vertices 

describing  fa 

ces.  Each  lis 

has  a flag  'if lag' 

as 

its  first  entry.  A 

value 

of  '1' 

indicates  that 

the  list  describes 

the 

edges  of  a face,  a 

value 

of  '0' 

indicates  that 

the  list  describes  edges  of  an  antiface.  The  next  entry  in  the  list  is 
' nv'  which  gives  the  number  of  vertices  in  the  list.  Following  ' nv' 
there  are  the  indices  ' i ’ of  the  vertices.  The  first  index  is  repeated 
as  the  last  one.  Note,  that  lists  describing  antifaces  need  not  follow 
the  list  describing  the  corresponding  face  but  may  be  placed  anywhere  in 


the  list  of  faces. 


’he  format  of  this 


ist  is: 


N 

13 

iflag1 

nv^  abed 

2013 

iflag2 

nv0  e b a 

2013 

iflagN  nv^  f g h 2013 

Note,  that  for  a pair  of 
must  he  one  and  only  one  pair  ' j i ' . 


(number  of  faces 
(first  face 
(second  face 

(last  face 

indices  of  vertices  'i  j’,  there 


5 . 3 Format  of  Vantage  Point  Specifications 

The  type  of  drawing  to  be  obtained  from  a vantage  point  and  the 
coordinates  of  the  vantage  point  given  in  the  object  coordinate  system 
are  described  by  ITYPE,  X0,  Y0,  20  on  a single  card  image  which  is  read 


format  ( IX , 

Ih,  5X,  3F10. 

3)  where 

ITYPE  = +1 

orthographic 

projection 

; all  lines  drawn 

+2 

orthographic 

projection 

; hidden  lines  dashed 

+3 

orthographic 

projection 

; only  visible  lines  drawn 

-1 

perspective 

projection-, 

all  lines  drawn 

-2 

perspective 

projection ; 

hidden  lines  dashed 

-3 

perspective 

projection; 

only  visible  lines  drawn 

X0  , Y0  , Z0 

are  the  x , y and  z coordi 

nates  of  the  vantage  point 

in  the  object 

coordinate 

system. 

The  project 

ions  and  the  coordinate  s, 

ystems  used  by  the  oregran 

described  in  Appendix  A. 
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5 . U An  Example  of  a Data  File 

A sample  data  file  of  a scene  is  listed  in  Fig.  5.1.  The  scene 
consists  of  three  objects:  a cube,  a pyramid  and  an  octahedron. 
Following  the  last  object,  two  vantage  points  are  specified.  The 
views  obtained  from  these  vantage  points  are  drawn  in  Fig.  5.2-3. 
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6.  THE  CDC  6600  BATCH  IMPLEMENTATION 

In  this  section  we  describe  system  programs  and  data  structure 
modifications  which  apply  only  to  the  CDC  6600  version  of  the  program. 

A complete  guide  explaining  how  to  run  PDRAW  on  a CDC  6600  and  how  to 
modify  it  for  use  on  similar  computers  is  being  issued  as  a User 
Bulletin  [ 5 ] ■ 

6.1  Description  of  CDC  6600  System  Programs 

In  this  section  three  CDC  6600  system  programs  required  by  PDRAW 
are  described. 

6.1.1  Subroutine  PLOTS 

Subroutine  PLOTS(LIMIT,  TEXT)  is  used  to  initialise  the 
plotting  output.  It  is  called  before  any  other  plotting  subroutine. 

LIMIT  is  the  total  length  of  a plot  in  integer  inches.  TEXT  is  a Hollerith 
string  of  up  to  ^0  characters  identifying  the  user.  PLOTS  is  a part  of 
the  CDC  66OO  SCOOP  plotting  package. 

6.1.2  Subroutine  PLOT 

Subroutine  FL0T(X,  Y,  IPEN)  is  used  to  control  all  pen  move- 
ments, create  an  origin  and  terminate  the  use  of  the  plotting  package. 

X and  Y are  coordinates  (in  inches)  to  which  the  pen  is  moved  from  its 
present  position  relative  to  the  current  origin.  IPEN  car.  have  these 
values : 

1 draw  a dashed  line 

2 draw  a solid  line 

3 move  the  pen  in  the  up  position 
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-1,  -2,  -3  the  pen  is  moved  as  for  1,  2 and  3 

respectively.  In  addition  the  final 
position  of  the  pen  (X,  Y)  becomes  the 
new  origin  (0,  0) 

999  terminates  the  plotting  output. 

PLOT  is  a part  of  the  SCOOP  plotting  package. 

6.1.3  Subroutine  SECOND 

Subroutine  SECOND  (TIME)  returns  the  elapsed  CPU  time  in 
seconds.  SECOND  is  a part  of  the  CDC  FTN  Fortran  Library. 

6.2  Notes  to  the  CDC  6cC0  Data  Structure 

The  following  notes  apply  to  the  CDC  6600  implementation  of  the 
data  structure: 

1.  The  arrays  in  the  data  structure  common  blocks  are  dimensioned 
as  follows : 


MAXM. 

20 

(the 

maximum  number  of 

cbj  ects ) 

MAXF . 

200 

(the 

maximum  number  of 

faces ) 

MAXE . 

800 

(the 

maximum  number  of 

edges  *2) 

MAXV. 

200 

(the 

maximum  number  of 

vertices  +1 

MAXS . 

Uo 

(the 

maximum  size  of  edge  intersection  stack 

This 

requires 

the 

following  CM  fields 

required  to  load 

661 ooQ 

required  to  run 

5^200^ 

program  sizes  are: 

FDRAW  programs 

10000- 

Common  31ocks 

2U500S 

system  programs 

20000,, 

- 3U  - 

The  dimensions  of  arrays  in  common  blocks  car.  be  easily  charged 
by  using  the  text  editor  UPDATE  and  a common  deck  containing  all 
common  block  statements. 

2.  The  following  numbers  were  found  to  be  good  computational 
tolerances  EPS1 , EPS2  and  EPS3:  0.1,  0.01,  l.On-iO. 

6. 3 Examples 

The  program  was  tested  with  scenes  of  one,  two  and  three  aircraft 
[Pig.  6.1-2],  Each  aircraft  consists  of  about  3^0  vertices,  330  faces  and 
660  edges . 

The  following  table  gives  memory  requirements  and  execution  times.  The 
packed  CDC  version  packs  integer  data  into  60-bit  memory  words  as  described 


Scene 

Unpacked  memory 

Packed  memory 

_ iT.e  " see 

3 aircraft 

210k 

115k 

11-25 

2 aircraft 

150k 

7^k 

6-15 

1 aircraft 

110k 

53k 

Example  [Pig.  5.2-3] 

33k 

33k 

0.01-0.05 

7.  THE  ADAGE  AGT-30  INTERACTIVE  IMPLEMENTATION 


7.1  Review  of  the  Implementation 

It  was  decided  to  implement  the  algorithm  on  an  Adage  AGT-30 
interactive  graphics  computer.  The  hidden-line  scenes  are  displayed 
on  a refresh  CRT  terminal.  Hard  copies  of  the  drawing  may  be  made  on 
a dot  matrix  plotter  or  on  an  incremental  plotter.  The  user  can 
interact  with  the  program  through  the  use  of  interactive  graphics  I/O 
devices.  The  Adage  computer  also  has  a good  Fortran  IV  compiler 
which  produces  efficient  code,  and  an  easy-to-use  display  text  editor. 

A major  problem  encountered  in  writing  the  program  was  the  limited 
amount  of  available  main  core  memory.  The  Adage  has  l6k  of  30-bit  words 
of  core  memory,  of  which  about  12k  is  available  to  the  user.  It  was 
decided  to  make  the  program  entirely  core  resident  for  minimum  execu- 
tion time.  This  is  important  to  facilitate  interactive  graphic  design. 

The  entire  data  structure  of  the  program  is  stored  in  several 
common  blocks  and  divided  into  two  parts:  vantage-point-independent 
data  and  vantage-point-dependent  data.  The  vantage-point-dependent 
data  is  stored  in  three  buffers  of  the  Adage  operating  system  AMOS. 2. 5. 
The  three  buffers  are  available  to  the  user  when  the  computer  does  no*: 
perform  any  disk,  tape  or  printer /'plotter  I ’0  operations.  This  save! 
about  2k  of  user  core  memory.  Another  saving  of  about  1.2k  of  memory 
was  made  by  simplifying  the  graphics  processor  3RAFX  into  a version 
which  displays  only  two-dimensional  vector  lists  and  ices  not  transom 
images  except  for  scaling. 
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The  vantage-point  coordinates  are  computed  from  sampled  values  of 
either  the  joystick  or  three  variable-control  dials.  An  object  des-. 
cription,  called  a model,  is  permanently  stored  in  a disk  file  and 
entered  into  the  program  by  a keyboard  command.  This  allows  the  user 
to  specify  interactively  a scene  and  a vantage  point.  A complete 
guide  to  running  PDRAW  on  an  Adage  computer  is  available  in  the  form 
of  a User  Bulletin  [6]. 

| 

T . 2 Execution  Time 

The  execution  time  of  hidden-line  removal  increases  with  the 
complexity  of  the  scene.  For  a simple  scene,  consisting  of  a single 
convex  polyhedron,  the  execution  time,  including  the  time  spent  by 
refreshing  the  CRT,  is  about  50  ms.  A model  of  a cube  can  "follow"  the 
movements  of  the  joystick  in  "real"  time.  For  more  complicated  scenes 
the  execution  time  increases  to  several  seconds.  In  order  to  avoid  a 
blank  CRT  screen,  during  hidden-line  computations  the  output  buffers 
for  display  lists  were  partitioned  into  two  parts.  While  lists  in  one 
part  are  being  displayed,  the  program  computes  new  display  lists  into 
the  other  part  of  the  display  buffers.  When  the  program  finishes  the 
removal  calculations,  the  display  processor  is  switched  to  display  "he 
new  lists. 

7. 3 Description  of  Interactive  PDRAW  Programs 
7.3.1  Program  PDRAW 

The  main  program  in  the  interactive  version  is  quite  differ- 
ent from  the  main  program  in  the  batch  version.  This  pr  gram  samples 
the  function  switches  and  branches  to  various  sections  upon  the 
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depression  of  one  of  them.  It  also  samples  the  variable  control  dials 
and  the  joystick  and  computes  the  coordinates  of  the  vantage  point  from 
their  values . PDRAW  includes  two  image  subprograms , PICTURE1  and 
PICTURE2,  which  alternate  in  displaying  two  sets  of  image  lists.  PDRAW 
starts  the  modified  display  processor  XGRAEX  by  a DISPLAY  statement  and 
stops  it  by  a call  to  NHALT  which  is  an  entry  point  into  the  display 
processor . 

7.3.2  Function  ALINE 

ALINE  is  a real  junction  which  returns  as  its  value  the 
y-axis  intersection  of  an  edge.  It  also  returns  the  slope  of  the  edge. 
ALINE  is  called  twice  by  ITEST  to  help  in  determining  the  intersection 
of  two  edges.  In  the  batch  CDC  6600  version,  this  function  is 
inserted  in-line  in  function  ITEST;  in  this  version  of  the  program  it 
is  made  into  an  external  function  in  order  to  save  core  memory. 

7.3.3  Function  CROSS 

This  is  a real  function  which  returns  the  magnitude  of  the 
cross-product  of  two  three-dimensional  vectors.  It  is  called  from  1 
locations  in  function  INTERS.  The  reasons  for  making  CROSS  into  an 
external  function  are  the  same  as  those  for  function  ALINE. 

7.3.1*  Subroutine  VECTOR 

This  subroutine  is  called  in  the  interactive  version  where 
calls  to  subroutine  PLOT  to  draw  a vector  are  made  in  the  batch  version. 
VECTOR  converts  the  coordinates  of  an  edge  or  a segment  of  ar.  edge  ir.tr 
a 2-D  display  format  and  appends  them  to  a display  list  -f  either 
visible  or  invisible  vectors. 
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7-1*  Description  of  Interactive  System  Programs 

In  this  section  three  interactive  system  programs  required  to  run 
PDRAW  on  an  ADAGE  A0T-30  computer  are  described. 

7.4.1  Program  XGRAFX 

This  program  is  a simplified  version  of  the  graphical  pro- 
cessor GRAFX.  XGRAFX  displays  only  lists  of  two-dimensional  vectors  and, 
except  for  scaling,  does  not  perform  any  image  transfonnations . These 
simplifications  save  about  1.2k  of  user  main  memory.  The  display  pro- 
cessor has  an  entry  point  MALT  which  stops  the  display  of  an  imago. 

It  is  called  before  a call  to  INPUTM  to  stop  the  display  while  INPUTM  is 
reading  from  a disk. 

7.4.2  Program  XPUIICHX 

This  program  punches  the  currently  displayed  image  in  a 
special  format  on  paper  tape  via  the  fast  paper  tape  punch.  The  paper 
tape  then  may  be  used  to  plot  the  image  on  an  incremental  digital 
plotter  connected  to  a PDP-8  computer.  The  PDP-8  plotting  program  is 
called  PLOT. 

7.4.3  Function  TIMER 

This  function  measures  elapsed  execution  time.  There  are 
two  entry  points  in  TIMER  called  TIMERON  and  TIMEROFF. 

7. 5 The  Data  Structure 

7.5.1  Common  Block  OUTPUT 

This  common  block,  used  only  in  the  Adage  interactive 
version,  contains  lists  cf  ail  the  visible  and  hidden  vectors  generate! 
by  the  program.  These  lists  are  in  a format  suitable  for  iiroo4:  iisplay 


P 
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on  the  Adage  CRT  by  the  modified  display  processor  XGRAFX. 

IVIPNT  is  a pointer  to  the  next  available  location 
in  a list  of  visible  vectors. 

IHIPNT  is  a pointer  to  the  next  available  location 
in  a list  of  hidden  vectors. 

These  two  pointers  are  initialized  by  program  PDRAW 
and  incremented  in  subroutine  VECTOR. 

IVISIBLE  (MAXVI ) is  an  array  containing  two  lists  of 
visible  vectors. 

IHIDDEN(MAXKI ) is  an  array  containing  two  lists  of 
hidden  vectors. 

Arrays  IVISIBLE  and  IHIDDEN  are  each  partitioned  into 
two  parts;  while  lists  in  one  part  of  the  arra2/s 
are  being  displayed,  the  program  computes  new 
lists  into  the  other  part  of  the  arrays  and  then 
switches  the  display  to  these  new  lists.  The 
visible  vectors  are  displayed  as  solid  lines,  the 
hidden  vectors  can  be  displayed  as  dashed  lines. 

'.5.2  Notes  to  the  Adaga  Tata  Structure 

The  following  notes  apply  to  the  Adage  AGT-30  implementation: 

1.  The  arrays  described  kc^this  pait  which  contain  vantage- 
point  dependent  information  are  put  into  three  common 
blocks  MONBUFFER 1 , M0MBUFFSR2 , and  MONBUFFER 3 . These 
common  blocks  are  stored  in  three  overlay  buffers  of 
the  Adage  operating  system.  In  this  way  they  share 
space  with  I ’0  drivers  and  save  ab  ".it  Ik  of  user  cere 


memory . 


_ 
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2. 


5. 


6. 


The  arrays  in  the  data  structure  common  blocks 


are  dimensioned  as  follows: 

MAXM 10  (maximum  number  of  objects) 

MAXF 50  (maximum  number  of  faces) 

MAXE  200  (maximum  number  of  edges  * 2) 

KAXV  50  (maximum  number  of  vertices  + l) 

MAXS  20  (maximum  size  of  edge  inter- 

section stack) 

MAXVI  ....  300  (maximum  size  of  visible  vector 
display  lists) 

MAXHI  ....  200  (maximum  size  of  hidden  vector 
display  lists) 


The  subroutine  TEAMS  scales  down  all  vertices  in  the 
picture  plane  (array  V2D)  to  (+1,-1)  square  in 
order  to  make  them  suitable  for  display  on  the  Adage 
CRT.  The  point  at  "infinity"  stored  in  V2D  (MAXV ,1 
and  V2D(MAXV,2)  is  set  to  (l.l,  l.l). 

The  computational  tolerances  in  EFSi,  EPS2  and  EPS; 
ir.  soBBBOn  block  ? ARAMS  are  set  to  0.1,  O.'l  and 
0.000005,  respectively. 

SCALE  in  common  block  PARA.M3  is  set  to  1.0  but  not 
used  in  this  version. 

DASHED  in  common  block  P ARAMS  is  used  only  in  the 


main  program. 


8.  CONCLUSION 


A program  implementing  the  Loutrel  hidden-line  algorithm  has  been 
written  and  tested  on  a CDC  6600  in  a batch  version  and  on  an  Adage 
AGT-30  in  an  interactive  version.  The  program  was  extensively  tested 
on  both  computers  which,  however,  does  not  guarantee  the  absence  of 
any  remaining  bugs. 

Some  of  the  improvements  and  extensions  which  may  be  made  to  the 
program  are  the  following. 

In  an  implementation  of  the  interactive  version  of  the  program  on 
a computer  with  a larger  main  memory  further  improvements  may  be  made 
to  the  algorithm  as  suggested  in  [U],  The  problem  considered  is:  given 
a scene  and  a vantage  point  which  moves  around  in  on  a trajectory,  it 
is  not  necessary  to  apply  the  entire  algorithm  to  each  frame,  but  rather 
we  would  take  advantage  of  the  similarities  between  successive  frames. 

We  note,  that  the  scenes  in  Fig.  3.1  and  Fig.  8.2  are  topologically 
equivalent.  That  is,  each  face  has  the  same  visibility  and  each  edge 
has  the  same  intersections  with  ether  edges  in  both  scenes.  The  two 
scenes,  however,  are  not  equivalent  geometrically.  To  obtain  the  scene 
in  Fig.  3.2  we  need  to  recompute  the  picture  plane  coordinates  of  all 
irertices  an ! ;dge  intersections  from  Fig*  3.1.  The  view  of  the  scene 
in  Fig.  8.3  is  not  topologically  juiva]  nt  • the  vi  ws  t s ne 
in  Fig.  .1  Fig.  . : ve  re-apply  th<  i jorittan. 

method  would  significant ly  reduce  the  execution  * imo  f or  m st  fram-s 
and  should  be  used  in  an  interactive  version  where  the  -"cecuti  n 
time  is  considered  to  be  import 


anl  ■ 
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Other  improvements  may  include  the  shading  of  surfaces  and  drawing 
the  scenes  on  gray-level  display  devices.  The  program  could  he 
extended  to  allow  a vantage  point  to  he  "inside"  a scene,  that  is 
inside  an  object  or  between  objects. 

In  other  implementations  of  this  program,  subroutines  SETUP,  INPUTM 
and  PRECOM  may  be  made  into  overlays,  which  would  save  space  for  the 
data  structure. 


L 
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Appendix  A 

THREE-DIMENSIONAL  TRANSFORMATION  AID  PROJECTIONS 

We  describe  here  the  transformation  from  a three-dimensional 
coordinate  system  to  a two-dimensional  picture  plane  coordinate 
system  as  well  as  the  two  projections  used  by  the  program. 

We  define  the  following  rectangular  coordinate  systems 
[Fig.  A. l] : 

1.  The  object  coordinate  system  (x,y,z)  in  which  the  vertices 

of  all  objects  as  well  as  all  vantage  points  are  defined.  The  vertex 

coordinates  are  entered  in  the  program  from  input  data  files.  The 

coordinates  of  a vantage  coint  specified  by  Q(x  ,v  ,2  ) are  ccmruted 

000 

from  sampled  values  of  either  the  joystick  or  variable-control  dials. 

The  object  coordinate  system  is  a right-handed  system. 

2.  The  vantage-point  coordinate  system  (x',y',z')  is  a left- 
handed  system  with  origin  at  the  vantage  point.  The  z'  axis  goes  from 
Q(xQ,yo,zo)  through  the  origin  of  the  (x,y,z)  system,  and  the  x'  axis  is 
parallel  to  the  x-y  plane  of  the  (x,y,z)  object  coordinate  system. 

3.  The  picture  plane  coordinate  system  (x",y"  is  a two- 
dimensional  coordinate  system  whose  axes  lie  in  a plane  that  is  perpen- 
dicular to  the  z'  axis  and  goes  through  the  origin  of  the  x,y,z 
system.  The  x"  axis  lies  in  the  x-y  plane  of  the  . x,y,z  system. 

In  the  program,  the  coordinates  of  the  object  system  x,y,z  . re- 
stored in  array  VERTEX . The  coordinates  f t le  -•  Lai  sy  • sc 

(x',y',z')  are  stored  in  array  VTRANS.  They  are  saved  only  in  crier 
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to  speed  up  computations  during  hidden-line  removal  in 
Coordinates  of  the  picture-plane  coordinate  system  (x" 
in  array  V2D . They  are  used  for  display  as  well  as  in 
intersections . 

The  transformation  from  the  object  system  (x,y,z) 
point  system  (x',y',z')  consists  of  the  following  four 
1.  Map  Q(x^,yo,ZQ)  into  the  origin  of  the  (x' ,y' 


T 


1 


1 

0 

0 

0 


0 

1 

0 

0 


0 

0 

1 

0 


-X 

o 


Xj 

o 
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[Fig 


2.  Rearrange  the 

8.2]: 


axes 


in  to  a left-handed  coordir 


0 

0 

0 


0 


0 


0 

1 

0 

0 


0 

0 

0 


3. 


Rotate  about  the  y'  axis  by  angle  a 
2 2 

d = x"  + y1"  cos  a = y ,’d 

0 0 o 


[Fig.  A." 
sin  cx  = 


function  INTERS. 
,y")  are  stored 
computing  edge 

to  the  vantage 
steps : 

,z ' ) system: 


ate  system 


i . 

J • 

x /d 
o 


U.  Rotate  about  the  x'  axis  by  angle  3 [Fig-  A.i]: 


The  four  matrices  can  be  multiplied  into  a single  transforms 
ion  matrix:  T = T,  T3  T. 


We  obtain: 
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The  transformation  from  the  vantage-point  coordinate  system 


(x',y',z')  to  the  picture-plane  coordinate  system  (x",y")  depends  on 


the  projection  system  that  is  selected: 


1.  Perspective  projection: 


In  this  projection  we  have  to  transform  using  these 


two  equations  [Fig.  A. 5]: 


x"  = D x'/'7' 


y"  = Dy' /z' 


2.  Orthographic  projection: 


In  this  projection  we  do  not  divide  by  the  depth  of  each  point. 


"he  equations  are  [Fig.  A. 6]: 


y"  = y 
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